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An  efficient  technique  of  microcomputer  mo¬ 
deling  is  presented  for  designing  linear  and 
nonlinear  circuits.  The  computer  model  is  valid 
for  pulse  forming  networks  (PIN's)  which  feed 
both  linear  and  nonlinear  loads.  These  PFN's  may 
resenble  classical  voltage  or  current  fed  one 
terminal-pair  Guillanin  networks,  but  are  in 
fact,  the  general  case  of  loaded  two  terminal- 
pair  networks  including  most  nonuniform  limped 
transmission  lines.  Before  building  systems 
using  PFN's  it  is  desirable  to  calculate  their 
response,  however,  closed  form  analytic  solutions 
are  impractical  and  other  computer  approaches 
such  as  SCEPTRE  and  SPICE  require  a  computer  with 
a  large  capacity.  An  algorithm  is  presented  that 
can  be  applied  to  pulsed  discharge  circuits  em¬ 
ploying  PFN's  and  solved  on  a  microcomputer  such 
as  an  Apple  or  HP-85.  Results  are  presented  for 
a  circuit  consisting  of  an  initially  charged  ca¬ 
pacitor  with  stray  inductance  in  series  with  a 
PFN  discharging  into  either  a  linear  or  a  non¬ 
linear  load. 
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Many  programs  have  been  developed  for  netr 
work  analysis  of  pulsed  circuits  such  as  SCEPTRE 
that  require  a  main-frame  computer  to  handle  the 
program  and  are  expensive  to  run  because  of  th^ 
time  needed  to  obtain  a  solution.  The  technique^ 
described  in  this  paper  has  the  advantage  of 
being  efficient  enough  to  be  programmed  on  a  mi¬ 
crocomputer.  The  microcomputer  is  more  readily 
available  and  much  less  expensive  to  operate  and 
produces  a  good  approximation  to  the  desired  wa¬ 
veform  in  a  short  solution  time. 
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where  tQ  =  initial  time,  and 

t„  =  last  time  the  network  is  to  be  solved 

r\ 

The  differential  equation  for  the  remaining 
loops  is  shown  in  Equation  2. 
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for  i=  2,3,4,.  ...n  (2) 

With  these  equations  and  knowing  all  the  ini¬ 
tial  values  of  current  and  voltage  at  time  t  a 
numerical  approximation  can  be  made  for  each  cur¬ 
rent  at  the  next  time  interval.  Eulers  method  ’ 
is  employed  which  uses  the  first  two  terms  of  a 
Taylor  series  expansion.  The  general  equation  for 
this  approximation  is  shown  in  Equation  3. 
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With  the  nw  values  of  current,  the  capacitor 
voltages  for  the  next  point  in  time  can  be  found 
by  numerical  integration.  The  general  equation  is 
shown  in  Equation  4. 
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for  i=  1,2,3, ... ,n  (4) 

This  methodgof  integration  differs  from  the 
rectangular  rule  in  that  it  approximates  the  area 
with  the  value  of  the  function  at  time  t^_^  multi¬ 
plied  by  the  increment  size. 


Given  the  n-section  LC  network  in  Figure  1, 
a  first  order  differential  equation  is  found  for 
each  current  by  suriming  the  voltage  drops  around 
each  loop  and  solving  for  the  derivative  of  the 
current.  The  differential  equation  for  the  first 
loop  is  given  in  Equation  1, 


It  is  assumed  that  the  voltage  across  the 
load,  v  +^,  is  defined  as  a  function  of  any  known 
value  of*  parameter  in  the  network  at  time  t^  or 
before,  for  instance; 
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These  are  sane  possible  forms  of  defining  the 
load  voltage. 

With  the  case  of  a  purely  resistive  load  the 
solution  is, 


n+1 


=  R„  x  i 


t=t,. 


Once  all  the  currents  and  voltages  are  found 
for  time,  t.,  the  process  is  repeated  for  time, 
t2  and  so  on. 

In  most  cases,  defining  loop  currents  (as 
opposed  to  branch  currents)  and  following  the 
paths  as  was  done  for  the  network  in  this  sec¬ 
tion,  gives  the  least  canplicated  and  shortest 
number  of  equations. 

QUOMIED  RESKNSE 

Taking  the  network  in  Figure  2  as  an  exam¬ 
ple,  the  differential  equations  for  current  i,, 
i 2 »  and  *3  acquired  by  sunning  the  voltage  drops 
around  the  loops  are  shown  in  Equations  5,  6,  and 
7,  respectively. 
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The  approximation  of  each  current  is  given 
in  Equations  8,  9  and  10. 


The  approximation  of  v.,  v„,  and  v,  is  given 
in  Equations  11,  12,  and  13,  respectively. 
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For  the  linear  case  of  =  R?  &(a  classical 
voltage  fed  Guillanin  type  A  network)3  the  voltage 
across  the  load,  v4  is  given  in  Equation  14. 


v4  t=tk  =  R2  x  ’l 


t=t. 


(14) 


This  algorithm  was  programed  on  an  HP-85  for 
parameter  values  designed  to  create  a  10  ns  pulse 
with  a  10  ohm  resistive  load  and  a  charging  vol¬ 
tage,  E  of  50  volts.  Other  than  the  charging 
voltage,  all  values  of  current  and  voltage  at  time 
t0  are  zero  and  the  time  incrsnent,  At  is  O.lus. 


To  verify  the  technique  an  experimental  model 
was  constructed  and  its  load  voltage  was  recorded. 
Figure  3  shows  the  canputer  and  experimental  vol¬ 
tage  waveforms  across  the  linear  load.  Error  due 
to  the  canputer  technique  is  of  interest,  however, 
it  was  difficult  to  distinguish  between  this  error 
and  errors  due  to  stray  inductance  and  capaci¬ 
tive  coupling. 


The  accuracy  of  this  technique  depends  on  the 
time  increment,  however  with  an  incronent  of  only 
l/100th  of  the  pulse  width  the  solution  takes  ap¬ 
proximately  44  seconds  and  the  percent  error  on 
the  average  of  the  amplitudes  of  the  first  five 
peaks  and  dips  of  the  waveform  in  Figure  3  was 
less  than  69s. 
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figure  3.  Computer  &  Experimentaiyollage 

Across  ih£  hi  near  Load  iM  ill 

In  order  to  further  demonstrate  the  capabi¬ 
lities  of  this  technique,  the  computer  waveform 
shown  in  Figure  3  is  canpared  to  the  waveform  ob¬ 
tained  by  analyzing  the  same  network  using  the 
SCU’TRE  program.  This  comparison  is  shown  in 
Figure  4. 
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The  root  mean  square  error  of  the  deviation 
between  SCEPTRE  results  and  results  obtained  by  the 
technique  described  here  is  ±  1.07  Volts. 
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FIGURE  4.  Comparison  Of  The  Technique  To  SCEPTRE 

The  output  waveforms  of  this  network  with  a 
nonlinear  load  consisting  of  a  100  ohm  resistor 
shunted  by  10  series-connected  1N914  diodes  were 
experimentally  recorded.  TTie  voltage  waveforms 
are  shown  in  Figure  5  and  the  current  waveforms 
are  shown  in  Figure  6 . 


FTCHJRE  IL.  Computer  A  Experimental  Voltage 
Across  Die  Nonlinear  Load 


FTfllRF.  fL.  Computer  A  Experimental  Qirrenl 
Through  Din  Nonlinear  Load 


Di screpenc ies  between  computer  and  experimen¬ 
tal  waveforms  can  be  attributed  to  the  model  of 
the  nonlinear  load  itself,  which  was  developed 
only  from  do  characteristics  of  the  load,  along 
with  the  errors  mentioned  in  the  linear  case.  On 
the  average,  the  percent  error  of  the  amplitudes 
of  the  first  five  peaks  and  dips  of  the  waveform 
in  Figure  5  is  less  than  495  and  the  solution  time 
is  approximately  51  seconds.  Figures  7,  8  and  9 
show  the  output  voltage  and  their  solution  time 
with  At  =  0.01  ns,  0.15  ns  and  0.2  ns  respecti¬ 
vely.  It  should  be  noted  that  with  a  time 


increment  of  0.01  ps  the  percent  error  is  less 
than  435  and  is  only  slightly  lower  than  that  of 
At  =  0.1  ps,  yet  the  solution  time  is  5  minutes  and 
45  seconds  longer.  If  At  becomes  too  large  so 
that  in  that  time  period  the  voltage  changes  a 
great  amount  the  computer  waveforms  will  oscil¬ 
late.  When  this  occurs  the  solution  time  actually 
is  longer  than  the  runs  with  smaller  time  incre¬ 
ments  because  it  takes  a  longer  time  to  plot  the 
oscillations.  Th i s  occurs  when  At  is  greater  than 
0.8  ps  in  the  linear  case  and  greater  than  0.15  ps 
in  the  nonlinear  case.  The  reason  this  value  of 
At  is  so  much  smaller  in  the  nonlinear  case  presen¬ 
ted  here,  is  because  at  low  currents  the  nonlinear 
load  acts  as  a  100  ohm  resistor  and  consequently 
the  voltage  changes  ten  times  more  for  each  value 
of  current  than  in  the  linear  10  fi  case. 


solution  time:  6  min.  36  s 
time  increment:  0.01  us 


E1QUKE  1 jl.  Computer  A  Exper  imental  Vol  tage 
Across  Hie  Nonlinear  Load 


solution  time:  36  s 
time  increment:  0.15  us 


FIGURE  JL.  Computer  A  Exper imentnl  Voltage 
Across  Hl£  Nonlinear  Load 


solution  time:  42  s 
time  increment:  0.2  us 

EIOJRE  iU  Computer,  4  Experimental  Voltage 
Across  Hie.  Nonlinear  Load 
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A  technique  has  been  presented  of  microcom¬ 
puter  modeling  which  can  be  used  in  the  design  of 
linear  and  nonlinear  circuits.  This  modeling 
technique  provides  a  good  approximation  of  the 
desired  waveforms.  Its  greatest  utility  is  in 
selecting  component  values  for  a  particular 
pulser  circuit  and  being  able  to  see  the  results 
before  the  circuit  is  actually  built.  The  advan¬ 
tage  of  this  technique  is  that  the  circuit 
modeling  can  be  run  on  a  microcomputer  found  in 
most  laboratories  or  universities,  the  solution 
times  are  reasonable  and  the  program  is  simple. 


BEEERQCTS 


1.  Bowers,  Sedore,  SCEPTRE;  A  Computer  Program 
fflt.  Circuit  and  Systems  Analysis,  Prentice-Hall, 
Inc.,  1971. 

2.  Sarjeant,  W.  J.  and  Hill,  R. ,  "Finite  Element 
Transient  Circuit  Analysis  Techniques,"  to  be 
published. 

3.  Huddleston,  John  V. ,  Nimerical  Methods  Lit 
Digital  Cbmputers,  c  1980,  pp  10.1,  10.2  and 
12.4. 


4.  Wylie,  Barrett,  Advanced  Engineering  Mathema¬ 
tics,  Fifth  Edition,  McGraw-Hill  Book  Co.,  1982, 

pg.  268. 


5.  Glasoe,  Lebacqz,  Pulse  Generators,  McGraw- 
Hill  Book  Co.,  1948,  Chapter  6. 


